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We calculate eigenvector statistics in an ensemble of non-Hermitian matrices describing open quan- 
tum systems [F. Haake et al., Z. Phys. B 88, 359 (1992)] in the limit of large matrix size. We show 
that ensemble-averaged eigenvector correlations corresponding to eigenvalues in the center of the 
support of the density of states in the complex plane are described by an expression recently derived 
for Ginibre's ensemble of random non-Hermitian matrices. 



The statistical properties of eigenvector overlaps may 
have an important bearing on time evolution and deter- 
mine the sensitivity to perturbations of systems governed 
by non-Hermitian operators or matrices. In such sys- 
tems it is thus important to know the statistical proper- 
ties of the (left and right) eigenvectors. Despite this fact 
little is known about eigenvector correlations in general 
ensembles of non-Hermitian random matrices. In Refs. 
jl],D, eigenvector statistics were calculated for Ginibre's 
ensemble of random non-Hermitian matrices where each 
matrix element is an independent, identically distributed 
Gaussian complex random variable. The question arises 
to which extent these results are relevant for other en- 
sembles of non-Hermitian random matrices. 
In the following we determine eigenvector statistics for an 
ensemble of non-Hermitian N x N matrices of the form 

M 

J = H + i 7 ^> Q (^ a ) t . (1) 

a=l 

Here H is a N x N Hermitian random matrix with com- 
plex, Gaussian distributed matrix elements Hki with zero 
mean and variance (\Hki\ 2 ) = Ski N" 1 . The V a are vec- 
tors with complex random Gaussian entries V k a with zero 

mean and variance (V/j*Vj) = SkiS a bN~ l . The eigenval- 
ues, A Q , of J are distributed in the complex plane. The 
ensemble (|l|), with 7 < 0, has been used to model the 
statistical properties of resonances arising in the case of 
resonance scattering in open quantum systems [^J^]; the 
position of the resonances is modelled by the real part 
A^ of the eigenvalues of (|l]) and the width by the imag- 
inary part A^ . The statistics of eigenvectors for such an 
ensemble was found to be of considerable importance for 
describing the properties of random lasing media, see || . 
The ensemble-averaged density of states d(z) = 
(N^ 1 J2 a S(z—X a )) for the ensemble ([!]) has been worked 
out using a number of different techniques, namely the 
replica trick |3| , the non- linear sigma model approach |^] 
and using the self-consistent Born approximation (tJ . 
Very recently, in Ref. M , all n-point spectral correlation 
functions for a variant of the ensemble (0) were deter- 
mined. It is given by 

J = H + if (2) 



where H is defined as above, and T is a fixed, N x N 
diagonal matrix with M non-zero diagonal matrix ele- 
ments 7. For the ensemble @j it was shown, in partic- 
ular, that the spectral two-point function ^2(^1,22) = 
(A^ 1 Xi tt ^^i - A Q )c5(z 2 - (and all higher corre- 
lation functions) are - after suitable rescaling and suf- 
ficiently far away from the boundary of the support of 
the spectrum - identical to those derived for Ginibre's 
ensemble (see Eqs. (15.1.31) and (15.1.37) of Ref. @). 
One thus expects that spectral n-point correlations of 
the ensemble (0) are locally similar to those in Ginibre's 
ensemble and thus universal. Furthermore it has been 
argued that under very general circumstances the fluctu- 
ations of the ensembles ([!]) and (||) are identical |J . 
Below, we explore to which extent the statistical prop- 
erties of eigenvectors in the ensembles (|l|) and (|^) are 
universal and the remainder of this paper is organized 
as follows. After defining the eigenvector correlator to 
be calculated, we briefly discuss the method used: the 
self-consistent Born approximation. We then derive an 
expression for eigenvector correlations and compare it 
to results of previous calculations for Ginibre's ensemble 
. Finally, we show results of numerical simulations, 
compare them to our analytical results, and discuss the 
applicability of our analytical method. 
The eigenvalues of ([!]) are non-degenerate with probabil- 
ity one, and in this case the left and right eigenvectors, 
\L a ) and \R a ), 

J \R a ) = A Q \R a ) , (3) 
(L a \ J = (L a \ X a 

form two complete, biorthogonal sets, and can be nor- 
malised so that 

(L a \R fj ) =6 af3 . (4) 

We indicate Hermitian conjugates of vectors in the usual 
way, so that, for example, \L a ) satisfies \L a ) — 
A Q \L a ). We investigate the eigenvector correlators |j],|2j 

0(z) = (^J2°^ S ( z - X «))> ( 5 ) 

a 

0(21,22) = {ij2 °<*e 5 ^ Zl ~ Aq ) 6 ^ - ( 6 ) 
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where O a f3 = {L a \Lp) {Rp\R a ) . These quantities may 
be extracted from 



D(zi,z 2 ) = {^'^O al 3 5(z 1 - X a )5(z 2 - \ 



(7) 



a,/3 



which may be written as D(z%, Z2) = 0(z\) S(zi — Z2) + 
0(zi,z 2 ). An expression for the diagonal part O(zi) for 
the ensemble ([l]) was derived in flD] ]. 



(G) = G + G £(G) 



(10) 



where Go = H^ 1 and S is a self-energy. A graphical 
representation of the self-energy (valid for M, N large 
and M/N = m = const.) is given in Fig. [j] The dia- 
grammatic rules are analogous to those described in [||. 
Differences are briefly explained in Fig. |^. Eq. (10) is 
solved for (G) in the limit of 77 — > 0. Expressions for the 
averaged Green function are given in . 
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FIG. 1. Shows a diagrammatic representation of the 
self-energy X. For the diagrammatic rules see J2) and Fig. 
^| below. 

Self- consistent Born approximation. We calculate eigen- 
vector correlators in terms of averages of products of 
Green functions, using an approximation scheme, namely 
an expansion in powers of N . The method used in 
Refs. , yielding exact results for Ginibre's ensemble, 
is not readily generalizable. We use the approach de- 
veloped in [[7| p^JT^ |: since the Green functions are non- 
analytic in the lower (upper) complex half-plane, a Her- 
mitian 2N x 2N matrix H = Hq + Hi is introduced 
T2l 



H u = 
with n > 0, A 



H 1 — j t 

-r) I I A f 

J and with inverse 



A 



G = 



G\\ G12 
G21 G22 



(8) 



(9) 



The resolvents are obtained by taking rj — > 0. In this 
limit, G21 = (z — J) -1 and G12 = (z — J^)^ 1 - Expand- 
ing the Green function G as a power series in Hi, its 
ensemble average (G) can be written as 




(c) 



FIG. 2. (a) Diagrammatic representation of the variance 
of the matrix H in t\m . (b) Representation of the second term 
of J [Eq. (^]. The random variables Vj? are denoted by wavy 
lines, (c) Shows the self-consistent equations determining the 

vertex . There are three more such equations determining 
the three remaining vertices occuring in Fig. [j]. 



From (G21) one obtains the density of states d(z) in the 
usual fashion |5],[l2 13 1. In order to make connection with 
the results discussed in @], we specialize to the limit of 
small m. In this limit one obtains 



d(0, y) ~ const. + m/ (4-7ry 2 ) 



(11) 



for m/(g + 1) < 2y < m/(g — 1) and zero otherwise, 
compare Eq. (108) in Ref. ||. Here 2g = (7 + 1/7) and 
z — x + ly with x, y real. We will analyze eigenvector 
correlations in the center zq of the support of the density 
of states, zq — xq + iyo with Xq = and yo = m/(2g) 
where do = d(zo) ~ g 2 /{irm). 



FIG. 3. Vertex for calculating averages of products of 
Green functions in the limit of large M, N with M = mN. 



Eigenvector correlators. We make use of the relation 



1 d d 

D(zi,z 2 ) = -2^^—F(zi,z 2 ) 

7T Z OZl OZ2 



(12) 
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with 

F( Zl ,z 2 ) = (AT-iTr [{ Zl - J)- 1 (z 2 - Jt)-i] ^ . (13) 

An expression for this average may be derived as de- 
scribed in J2|; cf. this reference for a diagrammatic rep- 
resentation of F(z\,Z2). The only difference between 
the case of interest here and the one discussed in Ref. 
Q is that the vertex must be replaced by that shown 
in Fig. |^. The corresponding expression for 7^(21,22) 
is valid in the limit of M, N large, M = mN and for 
\z 1 ~ z 2 \ 2 >(Trd{z + ) N)- 1 with z + = (z 1 + z 2 )/2. 
Using Eq. we obtain a (rather lengthy) expression 
for 0(zi,z 2 )- To keep the formulae simple, we specialise 
further to the case where z\ and z 2 are in the vicinity of 
Zq and obtain to leading order in m and lowest order in 
|fa| 



0(z a ,z a + fa) ~ - I - — 7 



1 1 



2g 4 J 7T 2 \Sz\ 4 



(14) 



Corrections breaking rotational invariance are of higher 
order. Comparison with Eq. (8) of Ref. jjj shows that 
locally, near the center of the support of the density of 
states, the eigenvector correlations for the ensemble (|l|) 
are the same as those in Ginibre's ensemble, apart from 
an additional factor of m/g 2 — m 2 /(2g 4 ) which measures 
the strength of the non-Hermiticity in Eq. ([!]). Accord- 
ing to Eq. (|l4|) , eigenvector correlations are strongest for 
7 ~ 1, and vanish in the limit of 7 — > and 7 — > 00 
which correspond to symmetric and complex symmetric 
J, respectively. 

Given the expression (|l4]), we may estimate O(zo) up to 
a constant of order unity, in the way described in . We 
obtain, to lowest order in to, 



777 

O(z ) * N- 
T 



(15) 



which is consistent with the result derived in pTj ]. In 
Fig. [| we compare the expression (|l4|) to the full result 
for 0(21, z 2 ) - as obtained within the self-consistent Born 
approximation - and observe very good agreement for 
|fa| not too large. 

Eq. ( |l4| ) is valid to lowest order in |fa| _1 and provided 
\8z\ 2 > (ndoN)^ 1 . The behaviour of 0(20,20 + Sz) for 
smaller values of Sz may be understood as follows. As- 
suming 



0(20,20 + Sz) ~ (O a/3 ) R 2 {z 0l z + Sz) 



(16) 



for |fa| ~ |A Q — \p\ very small, one may estimate the 
two factors on the r.h.s. of this equation separately. 
First, if two eigenvalues A Q and A/3 of J are very close to 
each other, one may argue that the corresponding over- 
lap matrix element O a p scales as |A a — \p\ . This is 
seen by considering a 2 x 2 matrix J with arbitrary com- 
plex matrix elements Jki ■ Denoting its right eigenvectors 
by \R a ) — (l,p a ) t , the corresponding left eigenvectors, 



assuming Ai ^ X 2 and subject to the condition (||) of 
biorthogonality, are given by by {L\ | = (-Q2, l)/(ei - Q2) 
and (L 2 \ = (-Qi, 1)/(q2 - Qi)- Thus 



Ol 2 



OC 



(£l|i2><i?2|i2l) 

-|Ai-A 2 |- 2 . 



\Qi - Q2\ 



(17) 



For A Q very close to A/3 [namely on scales smaller than 
the mean level spacing (7rdo-/V) -1 / 2 ] this behaviour per- 
tains for arbitrary values of N . Second, the spectral 
two-point function scales as R 2 {z$, 20 + Sz) oc |fa| 2 when 
I fa I — * H||]. Thus one concludes that 0(zi,2 2 ) must 
converge to a constant as \Sz\ approaches zero. Since the 
crossover from Eq. jl4| ) to constant behaviour occurs at 
|fa| 2 ~ (ndoN)^ 1 , one estimates, to lowest order in m 



0(2o, 2: + fa) TO 2 



(d N)2 



fJ 



for |fa| 2 < (ird Ny 



(18) 



Eqs. (^4J) and (|18|) are consistent with the assumption 
that eigenvector correlations in the scattering ensemble 
(|l|) are universal in that 0(20, zq + Sz) is given, after suit- 
able rescaling and well within the support of the density 
of states, by the corresponding expression derived in Ref. 
|H for Ginibre's ensemble: defining ]li[] Sz = Sz^irdaN 
one may expect, to lowest order in m 



O(2o, 2p + fa) 

(d N)i 



(19) 



9 m 



[l-(l + |fa| 2 ) exp(-|fa| 2 )] 



This expression interpolates between (14) and (18|). 
It must be pointed out that Eqs. fll8]) and (|l|) cannot 
be valid for very small values of to = 0(1 fN), where 
the ensemble (||) deviates very little from the classical 
Gaussian unitary ensemble of random Hermitian matri- 
ces [|lO| . Spectral correlations in this situation have been 
analyzed in detail in Refs. |^|,^| where it was shown how 
the crossover from non-Hermitian to Hermitian ensem- 
bles may be characterized. 

Numerical results. We have verified the validity of Eq. 
( |l9| ) using numerical simulations of of the ensemble (Q), 
for to = 0.1,7 = 0.5 and N = 100,200,400,800 and 
1600. Fig. I shows 0(20, 2 + Sz)/{d N) 2 as a function 
of Sz. We observe that the numerical results converge 
to Eqs. ( p^| , p+| ). Convergence with increasing values of 
N is much faster for small values of |fa| than for large 
values of | Sz . In fig. [|, the scale of the a;-axis differs from 
that of fig. I (a) and differences between (|lj) and the full 
result - as obtained within the self-consistent Born ap- 
proximation - are not visible here. 

We have also performed simulations for the modified en- 
semble (0) . The results are very similar to those displayed 
in Fig. pi (not shown). 
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FIG. 4. Shows ReO(zo,Zo + Sx) as obtained within the 
self-consistent Born approximation (solid line) compared with 
the asymptotic expression @ (dashed line) (a) as a function 
of Sx for m — 0.1 and 7 = 0.5, (b) as a function of m for 
Sx — 0.01 and 7 = 0.5, (c) as a function of 7 for m = 0.1 and 
Sx = 0.01. 



Conclusions. In this paper we have calculated the eigen- 
vector correlator 0(zi,z 2 ) for the ensemble (||) using the 
self-consistent Born approximation, and for both ensem- 
bles (0) and (§) using numerical simulations. Our results 
imply that eigenvector correlations in these ensembles are 
locally given by a universal law - after suitable rescaling 
of the complex energies. One may thus expect that local 
eigenvector correlations in more general ensembles (such 
as ensembles of random Fokker-Planck operators 12|) 
may be described by the law derived in Q. It has been 
pointed out that such correlations may determine tran- 
sient features in the dynamics of such systems The 
results found here may be of direct relevance for quantum 
scattering systems P,|l6|]. 
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FIG. 5. Shows ReO(zo, zo + Sx)/(doN) 2 as a function 
of Sx = Sx y/ivdoN; for m = 0.1,7 = 0.5, N = 1600 (A), 
N = 800 (V), N = 400 (o), N = 200 (□) and N = 100 (O). 
Also shown is the analytical estimate according to Eq. (JIJ) 
(dashed line), and Eq. ( |l£| ) (solid line). 
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